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MOEH YDR- Snowme It Model 
For Augmenting Streamflow 
Prediction and Simulation 



L. A. Logan 



ABSTRACT 



A mathematical snowmelt model formulated on the basis 
of physical and empirical melt equations is developed in 
complement with a conceptual watershed hydrologic model 
(MOEHYDR) for daily streamflow prediction and simulation. 
Energy balance and aerodynamic principles provide the basis 
for the input/output heat energy requirement for the snow- 
melting processes, while the physical condition of the 
snowpack and its variations in depth and density are repre- 
sented through locally derived empirical equations. These 
equations are complemented in application, by more general 
literature-cited empirical equations that relate to other 
critical snowpack conditions and the environment. Four 
applicable data input options are formulated in recognition 
of the lack of availability of sufficient climatic data in 
most drainage basins. Computer implementation of the snow- 
melt model with one of the applicable data input options is 
accomplished as part of an overall MOEHYDR- hydrologic model 
simulation package* Test applications indicate that the 
snowmelt model is a satisfactory complement to the MOEHYDR- 
model in winter streamflow simulation in parts of southern 
Ontario. 



■-> ^ ..;-<■■-- 'm:^!^^'i>:^^^^F^->-'^^m-'"-'^:''''^m»:^%-: '^- '; "♦ -^ ■ ■■■" '^ 



MOEHYDR- Snowmelt Model for Augmanting 
Stxeamflow PredictiO'ii, and Simulatio'n 

L. A, Logan ^ 
INTRODUCTION 



Paxt of the activities of the Water Modelling Section 
O'f the Ontario^ Ministry of the Environment is involved with 
the development and implementation of hydrologic system 
models for use in the prediction or simulation of stream- 
flows in selected drainage basins in Ontario* In general ^ 
these models are employed as analytical tools and aid in 
solving specific water resources management problems. 

The sO'lution,s to many of the water management problems 
in hydrology require estimates of streamflow which are based 
on snowmelt computations. The simple temperature index 
method is often inadequate for computing reliable estimates 
O'f the excess snowmelt required for the synthesis of the 
streamflow. This is often and particularly true under 
southern Ontario winter conditions where, because of the 
usually shallow snowpack accumulations and large fluctua- 
tions in climatic conditions, the snowmelt rates within a 
given basin, may be influenced more directly and signifi- 
cantly by energy contributing variables other than tempera- 
ture. Therefore I the modelling approach considered here is 
based on the energy balance and aerodynamic principles. 
Consideration is also given to the known physical properties 
of the snowpack. This approach is similar to that employed 
by Amorocho and Espildora (1966) , with the introduction of 
some modifications to accommodate local snowmelting conditions 

A computer simulation package program for a total 
watershed hydrologic system model is developed in the Water 
Modelling Section of the Ontario Ministry of the Environment. 
The package program comprises, in addition to the snowmelt 
submodel, a number of other simulation submodels with the 
capability for both hourly and daily simulations. The 
objective of this paper is to present a general outline of 
the snowmelt submodel and sample results of its test 
applications. 

MOEHYDR - Watershed Hydrologic Model Package 

The model is developed and formulated to represent as 
best as possible, conceptualization of the various hydrologic 
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processes in the watershed that are realized from the 
precipitation input, changes in storages and evapotrans-- 
piration losses to runoff at the basin outlet. The model 
package comprises three major submodels* one for daily 
simulation, another for hourly simulation during specified 
storm periods and the other for simulating the effect of 
snow accumulation and snowmelt. All of the three submodels 
are linked and interfaced for sequential time discrete 
eimulations. The package program is completed with an 
automatic optimization routine and euiranary routines which 
comprise ten output options. Figure 1 outlines the general 
scheme of the model package, 

LTSM: Daily Simulation Submodel 

This submodel effects a time discrete, daily accounting 
of the watershed hydrologic processes, conceptually, using 
the principles of mass conservation and hydrologic contin- 
uity. The model comprises fifteen major accounting opera- 
tions (e.g. interception and depression storages, infiltra- 
tion, interflow, flow routing, etc,). Fourteen physically 
based parameters and eight coefficients and constants are 
employed. Four parameters and three coefficients are opti- 
mized automatically. The other parameters are estimated 
off-line and by trial and error optimization of the model 
under the guidance of sO'Und hydrO'logic principles and 
intuitions, 

STSM: Hourly Simulation Submodel 

This submodel performe a storm event based and hourly 
time discrete accounting of the watershed hydrologic pro- 
cesses, similar to that of LTSM, but on the finer time 
scale. This facilitates modelling of the small time basin 
responses. The model comprises twenty- two simulating and 
accounting operations. The smaller time increment facili- 
tates more refined, but more complex simulation algorithmi 
for the most sensitive and critical processes such as 
infiltration and overland flow. These algorithms are more 
realistic in the details of concepts of the mechanism of the 
hydrologic processes involved. Fourteen physically based 
parameters and ten coefficients and constants are utilized 
by this submodel. The parameters and coefficients of the 
modelled infiltration and overland flow processes are 
optimized automatically. The other supporting model para- 
meters and coefficients, optim,ized by the LTSM, are used as 
compiements in the overall sequential time discrete simulation 

SNMELT: Daily Snowmelt Simulation Submodel 

This submodel simulates for the winter mo^nths, the 
amount of snowfall lost to interception storage and evapora- 
tion or sublimation, increase in the snowpack depth and 
excess melt which goes as net moisture input for use by 
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FIGO'RE 1 . General Scheme of the 

MOEHYDR - HydrO'logic Model Computer 
Simulation Package . 



LTSM- The sno^wmelting algorithms utilize knowledge about 
the physics O'f the sno^'wmelt proces.ses and allow implicit 
consideration of the physical and thermal properties of the 
snowpack. The energy balance and aerodynamic approaches are 
used in relation to the eource of snowmelt energy supply. 
Six major simulation operations are executed, including 
simulation of snowmelt acco'rding to: no rain with negative 
and positive energy fluxes and rain-on- snow with negative 
and positive energy fluxes. The model utilizes a basic 
daily input of seven different meteorologic variables (soime 
as observed and others as computed or estimated inputs) . 
The submodel may operate in accordance with any one of four 
data input options provided to facilitate a range in the 
availability of local climatic data. 

Sets of locally derived empirically based equations 
involving snowpack density, depth and selected climatic 
variables, are employed in order to effect characteristic 

lO'Cal snowpack conditions in simulating the snowmelt 
processes. 

PATSEARHi Pattern Search Optimization Routine 

This routine utilizes the direct pattern search method 
of optimization (12) ^ . The technique comprises a pattern of 
moves, exploration and acceleration from given parameter 
base points. The procedure continues by updating the base 
points successively, with successful moves in relation to an 
improvement in minimizing the objective function with 
respect to the parameters. Four different forms of objec- 
tive functions may be used, namely sum of squared, weighted 
sum of squared, or absolute residuals, or the sum of predic- 
tions. The routine is controlled by the MAIN program which 
reads the variables to control the optimization procedure, 
optimizes the selected model parameters from given initial 
values and performs statistical analysis on the observed and 
simulated flows. The routine can optimize up to twenty 
parameters, and reproduces response surfaces for paired 
parameters. 

SUMMARY/ENDMON: Summary and Plot Routines 

These routines provide daily sunuiaries of the simula- 
tion outputs for all the significant model components and 
print out end of month summaries of storages and flows. In 
addition, eight alternative printout options (OPT 1, OPT 2, 

etc.) and three plot options (LINEP, OPT 7, OPT 8) are 
available. 



^ appended reference 
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THE; SNOWMELT (SNMELT) MODEL 



In the formulation of a physically based snowmelt 
model, appropriate consideration should be given to the 
supply of heat energy required for the snowmelting processes 
and also to the physical and thermal character ietics of the 

snowpack,. 

The basic equations that describe the heat exchange 
between the snowpack and its environment are well documented 
in several reports and texts (1,2,5,6,7,8). Other useful 
empirical equations are those that interrelate certain 
critical snowpack characteristic and the environment, such 
as density of snowfall in relation to air temperature 
(1|2,5) settlement or compaction of the pack (2,4) and 
density of the pack in relation to its depth and rainfall 
input (9, 10) . 

Physical Melt Equations; Melt Processes 

For a given snowpack of unit area, for any given time 
increment At, the heat equivalent of melt Bupplied from 
given heat sources, can be defined by way of an energy 
balance equation (5,7,8) i 

(1) HM.At ^ C-HC + HQ).At ... .... cal/cm'^ 

where ■ 

HM - the heat equivalent of melt (cal/cm^), 
HC = the heat capacity of the pack (cal/cm^ ) and 
HQ - the change in heat stored in the pack (cal/cm-); 
all for time increment At. 

The change in the pack's heat storage is effected by the 

external energy sources given by: 

(2) HQ.At = (HR + HA + HAC).At . ... cal/cm,^ 

where i 

HR = the net radiative (shortwave and longwave) 

energy flux (cal/cm^), 
HA = the net turbulent exchange (convection and 

condensation or eublimation) energy flux (cal/cm^) 

and 
HAC= the net advective (rain) and conductive (ground) 

energy fluxes (cal/cm^) ; all for time. At. 

The heat capacity HC, the heat required per unit area to 
raise the temperature of the pack to 0°C| is given by: 

( 3 ) HC . At = CS . (DENS . DEPS . TEMS ) . At ... cal/cm^ 



where i _ ^ , _ ^ 

CS = the specific heat of snow (%0.5 cal g ' ^'c ), 

DENS,, DEPS and TEMS axe the average density, depth 
and temperature, respectively, of the pack for 
time. At. 

The shortwave and longwave radiation are the itiajor sources 
of energy exchange with the melting snowpack. The net 
shortwave radiation, HRSN, that can be absorbed by a pack for 
any given time, At, may be expressed as; 

(4) HRSN. At = KT. (1-ALB) . (l-KC.CN.)CHRS.At cal/cm^ 

where i 

KT is the forest cover transmission coefficient, 
ALB is the snowpack albedo, 

EC is the cloud ceiling factor, 

CN is the degree of cloud cover, and 

CHRS is the clear sky direct insolation (cal/cm^); all 
for time, At,, 
Similarly, for any given time. At, the net longwave radiation. 
HRLN may be expressed by. 

(5) HRLN. At = (KF + (1-KF) .( 1-KCL. CN) ). SIGMA 

. ( EMA . TEMA "* - EMS . TEMS V) . A t ca 1/ cm_, ^ 

where i 

KF is the forest cover density factor, 

KCL is the cloud ceiling factor for longwave radiation, 

SIGMA is Stefan ' s constant ( . 826x10" ^ ° cal/cm^/min/OK'* ) , 

EMA and EMS are the eraissivities of air (0,759) and 

the snowpack (1.0), respectively, 

TEMA and TEMS are the temperature O'f air and the 

snowpack, in ^K, respectively for time. At. 

From a given height of Im above the pack, the convective 

heat exchange, HCONV, and the heat exchange by sublimation, 
HCOND, for any given time, At, may be expressed, respecti'- 
vely by: 

(6) HCONV. At = KCONV. ( (TEMA-TEMS) .VW) .At . . . cal/cm^ 

and 

(7) HCOND. At = LHS. (KCOND. (EA-ES) .VW) .At ... cal/cm^ 
where : 

KCONV and KCOND are transfer coefficients in accordance 

with the system of units i 

TEMA and TEMS are the air temperature at the fixed 

height and at the snowpack surface, respectively, 

for time, At F 

EA and ES are the vapour pressures of the air in 

mb, at the fixed height and at the snow surface respectively, 

at a given temperature, for time, At| 
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¥w is the average wind speed in km/h , at the fixed 
height for time, At, and 

LHS is the latent heat of sublimation or condensation 
(60 cal/g) . 

In the presence O'f rain-o^n-snow, an additional heat component 
of HQ is the advective heat, HP, which may be computed, for 
time. At, by: 

( 8 ) HP . A t = CW. DENW. PRN (TEMA-TEMS ) . At . . . cal/cm^ 

where i 

DENW is the density of the rainwater (~1.0 g/cm^), 
CW is the specific heat of water (1.0 aal/g/^C} and 
PRN is the net amount of rain-on- snow (cm) for 
time, At. 

The final component of HQ is the conductive heat across the 
ground-snow interface. The conductive heat HGD for time. 
At, may be expressed generally by: 

(9) HGD. At = KGD. (DTEM/DZ) .At ... . . . . . cal/cm^ 

where i 

KGD is the thermal conductivity of the soil and 
(DTEM/DZ) is the temperature gradient existing across 

the ground-snow interface for time, At, 

The potential meltyMT, in water depth equivalent for time. 
At, is given by (5,7,8): 

(10') MT. At = (HM + HC) / (TQ. LHFW) .At ... cm 

where : 

LHFW is the latent heat of fusion of pure ice at 0*^0 
(80 cal/g) and 

TQ is the thermal quality of the pack, at the start 
of the given time increment At, and is defined by: 

(11) TQ = LHFS/LHFW + (CS . TEMS) /LHFW 

where : 

LHFS is the latent heat of fusion of the snow. 

The snowmelting processes proceed according to the change 

in the heat storage of the pack. The value of HQ(>0 or <_0) 
and TQ(<1 or ^1) usually indicate the satisfying or non- 
satisfying of the pack's heat deficit, DM, given by: 

(12) DH.At - DENS. DEPS. (LHFS + CS. TEMS). At ...... cal/cm^ 

in which the change in TEMS is given by: 

DTEMS . At ^ ( HQ +LHF ) . A t/ ( CS . DENS , DEPS ) . A t . . . ^ G 
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where LHF is the latent heat of fusion releases by the freez- 
ing of the liquid water content^ WC, of the pack (if HQ'<0): 

(13)LHF . At = LHFW. DENS . DEPS . WC . At . . . cal/cm^ 

For TQ <1, WC may be expressed by: 

(14) WC.At = (1 - TQ) .At 

Excess melt is produced and is eventually drained from the 
pack when the melt exceeds the liquid-water holding capacity, 
WHC, of the pack 

(15) EMT.At = (WHC - WC) .DENS . DEPS .At .„., en! 

Empirical Equations- Snowpack Condition 

The empirical equations employed are related primarily to 

the cO'ndition of the snowpack associated with the pack's 
accumulatiO'n and depletio^n. Interception of snowfall by the 

vegetative canopy is approximate by: 

(16) PN.At = (1 - KV).P.At . ... cm 

where PN and P denote net and total snowfall respectively, 

and KV is an arbitrary constant for vegetation interception. 

The density of fresh snow is observed to be a function of 
air temperature. A suitable empirical equation for this 
relationship is given by (1, 2) : 

(17) DENEW At =-r°-°^ + (TEMA. At/100 )S for TEMA > 0°F 
(1/) DENEW. At - -| g_Q5 ^ gQj. ^j,j^ ^ qOj. 

Where : 

DENEW is the density (g/cm,^') of the fresh snow. 

As a result of the new snowfall the depth of the pack is 
increased, but, is also adjusted for compaction or settle- 
ment according toi 

(18) DEP.At = (WEQ + DENEW. DEPNEW) /DEN. At ... cm 

where : 

WEQ and DEP denote the initial water equivalent and 

adjusted snowpack depth, respectively, and DEN denotes 
the adjusted density of the pack computed by the 
empirical equation (1, 2) : 

(19) DEN. At = DENS(1 + KX. DENEW. DEPNEW) . At ... g/cm^ 

where i 

KX is a factor for change in density due tO' the compac- 
tion of the new snowfall and settlement. 
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similarly, the snowpack may be adjusted for eublimation (if 
dewpoint temperature <0°C) according to: 

(20) DEP.At - DEPS. (I-(DENW/DENS) . (HCOND/LHFW) ) .At .. cm 

where s 

DENW is the density of liquid water. 

HCOND is given by eq. (7) which provides the rnechanism 

for heat exchange with the pack due to eublimation. 

Empirical equations relating average snowpack density, depth 
and rainfall amount, are derived to complement the simulation 
of changes in the pack's density that are influenced by the 
prevailing local climatic condition. The relationships are 
derived for specified accumulation and melt periods, with 
and without rain-on-snow (9, 10): 



(21) DENS = A 4- B log (DEPiS) . g/cm'' 

(22) DENS = Al + Bl.log (DEPS') + B2.APR . g/cm' 

where i 



DENS, DEPS and APR denote the average density and 

depth of the pack and average rainfall input reBpectively 

Finally, the water holding capacity, WHC, of the pack may be 
updated empirically by (2) : 

(23) WHC = 0.025, DENS + 0.03, DENS < 0.4 
0.2. DENS = 0.04 , 0.4 < DENS < 0.6 
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Simulation Algorithms- Computer ImplementatiO'n 

The following texts and flow charts outline, generally^ the 

simulation algorithms that are computerized. 

Auxiliary Input Procedures - Figure 1 shows the specifications 
for the simulation relative to basin segments of eleyation- 

zo'nesr implying some provisions for spatial variability in 
the snowpack accumulation and melt rates . 

The routine provides for simulation options depending on the 
availability of climatic data. Subroutines are executed for 
computing interpolated and/or weighted average climatic data 
input for the model. 



ENTER ELEVflTION -ZONE PMYSIO'GRAPHIIC PARflMiETIEna 
READ AREAJKI, PCAREAIK), FARiStK)'. ISCM (K 1', IMELV IK), 
K°lll, MM. 



lENITEiR ELEVATION- ZO'NIE DENSITY - BEPTH EMiPIRICAL 
CO'EFFICIENTS ft eOiNSTANTS ' ACCMLN. ft iMELT PERIO'DS 
iREAD AAP(K1,B.AP1K), AH'AP(K), B:RAP(Kli, BRAP HKl, AMIPtKl 
liMiP(K),A;RMPm),BRMP(K).SRMlPIII(K)', KH.MM 



ENTER ELEVATIOIN- ZONE TIME -AREA SNOWCOVER OEPLN. 
READ' SCDdK.Jil. K-l,Mlil. J»il,NNi 



ENTER iLEVATION-ZONE INITIAL SNOWIPACK CQ'NDITIO'N. 

READ DC!p«K,.J|.D€NlK„J).WEQlK,J),VillS(K.Jl„MlliST|K,J)., 

WCiK,J|i,TQl(K,JI, TEM(K,.J|,WHCt!K,Jl,K.||.liM,Jl.|!.NN 



SPECIF'T SIMIULATION OPTIO'N- 
O'PTION 1.2,3' 0mR4 



ENTER MlETEOlROLOSICAL O'ATA IFO'R OPTION BY GBOUPEiD 

VARIABLES - E.G. O'PTION "l" 

READ PR(l,.J),PStl.J|. I-I,M Jl.1ll,NIN TEMDliI .J), EA(I, J) 

P£T(l,J),HHSill,Jl,HBHI„Jll, ALB|l,Jl,VW|il,J) 

TEMACI.J). I-I/M J-II,NN 



ENTER ANiCIILLAiHY MELT COEFFICIENTS ii CONSTAftlTS 
READ BS, KX.CS, KCOND, KCONiV. MOD 



eALL"DISTPPT' 



CALL "DISTEMIP" 



CALL "DISTWIND" 



CALL "DlSTOMDTA" 



H® 



Figure 1: Auxiliary Input Procedures for the Snowmelt 
^~~^~~ Simulation Algorithm 



'4«1. 



-,. ,,.,.^,ll^,i^gp^,^. ^.^^^J^^^^^^JJP^ 



irm- 



11 



Interception Storage and Snowfall Density - Figure 2 shows 
the continuation of the simulation algorithm. Allowances are 
made for the reduction in precipitation input (rain or snow)^ 
for interception storage and for computing empirically the 
density of fresh fallen snow as a function of air temperature. 



j^i 



K = l 



AREAS AREAIK) 
PCAREA = PCAREAtKj 
PARIS = FA RIS(K) 
ISCM-ISCM(K) 
MKLV-MELVtK) 



HCOND = hCOND(K,J) 
HCONV = HCONVtK,J) 
HGD' HOD(K,J) 



ACAEIS^ 0.0 
ACAESP =0 
AET -0.0 



AAP = AAP(K) 
BAP = BAP(K) 
ARAP- ARAPCK) 
BRAP = BRAP(K) 
BRAPI =BRAPI(K) 
AMP - AMP{K) 
BMP ^ BMPIK) 
ARMP = ARMP(K) 
BRMP = BRMPtK) 
BRMPI=BRMPI(Kl 



i_ 



DEPMDEP(K,J) 
DENMDEN(K,J) 
WEQNWEQIK.J) 
WCI^WC(K,J> 
WHCI -WHCIK.J) 
TEMI = TEM(K,J) 
VIS = VIS(K,J) 
MIST^MIST^K.J) 
VISO =VIS 
SCD- SCOtK.d) 
TQI - TQ(K,J} 



-M2I 



VIS<ISCM 




HR = HRN(K,J) 
PR = PR(K.J) 
PS =PS(K,J) 
TEMA-TEMAtK.J) 
TEMD - TEMD(K,J] 
EA=EAtK,J) 
PET- PET(K,J) 
VW =VW(K,J) 
P = PR + PS 




IS = BS*PR 
PRN=PR-IS 
VIS = VIS + IS 
MIST = MIST + IS 




DEP0»O.O 
DEN 0=0.0 
WEQ 0-0.0 
W C = 0.0 




r 



VIS'iSCM 




DEN 0^0,05 




IS = BS«PS 
PSN = PS-IS 
VIS -VIS + 1! 



DEN0-f(TEMAJ 



Figure 2: Simulation Algorithm for Snowpack Parameter 

Initialization, Snowfall Density and Interception 
Storage, 
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Snowpack Accumulatioii - Figure 3 continues the simulation 
algorithm for snow accumulation and some reduction in depth 
by compaction due to the weight of the new snow. The heat 
exchange and the computations for evaporation or sublimation 
from the interception storage ^re also described. 

0^ — I 




DEP0 = PSN/DEN 
WEQ0 = DEN 0*DEP0 
WC = 0.0 
WHC -0.03 
TEM0 = TEMA 
TQ = 1.0 



DEP 2 -0.0 
DEN 2 = 0.0 
WEQ2 = 0.0 
WC 2 =0.0 
WHC 2 -0.0 
TEM2 «0.0 
TQ2- 1.0 



DEP I = DEP 
DENI = DEN 
WEQ 1 = WEQ 
WC I = WC 
WHC I = WHC 
TEM I = TEM 
TQ I = TO 



DEP -DEP 1 + DEP 

DEN = DEN I nil + KX « WEQ | 

DEPI=^DEP#DENI/p£f^ 

DENI =DEN 

WEQ I -DEN I* DEP I 

ivHCI = f (DENI) 

TQ I = l-5/g »(TEMI -32.0)m CS 



PEIS = PET« FARIS 



80 



1 



VIS = VIS -AEIS 
MIST - RMIS 
DViS = VIS -VIS 
ACAEIS s ACAEIS +AEIS 
AET = AET +AEIS 




<8) 



HQ=HR+ HCOND + HCONV 
MIS= HQ-t^FAKIS 

2.54*80 
MIS = MIST + MIS 



MIS *0.0 

PHE=ABS( HCOND) 
PSE= PHE» FAR IS 
2.54 *80.0 




AMIS ^ MIS 



AMIS -VIS 




AEIS =PSE 



1 



AEIS - VIS 



AEIS '^ PEIS 










RMIS = AMIS - AEIS 



AEIS- AMIS 



RMIS = 0.0 



L 



VIS = VIS - AEIS 



DVIS - VIS - VISO 
ACAEIS = ACAEIS + AEIS 
AET = AET + AEIS 



I 



kD 



Figure 3 : Simulation Algorithm for Snowpack Accumulation 
and Sublimation from Interception Storage. 
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Sublimation and Evaporation - Figure 4 outlines the algorithm 
for computing sublimation or evaporation from the snowpack, 
summation of the heat energy and the executions for the 
snowmelting processes for HQ, with or without rainfall on 
the snowpack. 
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Figure 4 : Simulation Algorithm for Sublimation/Evaporation 
from Snowpack. 
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Snowmelt - No Rain, Negative HQ - Figure 5 shows the simulation 



algorithm for "Snowme It-no- Rain Negative Energy Flux". 
Thus, at the start of At, if TQ <1, this implies that WC >0 and 
DH = 0. Therefore the HQ(-) then freezes all or portion of 
WC and the temperature of the pack is changed accordingly, 
depending on the amount of latent heat of fusion released by 
the freezing of WC. 

If TQ >1, it is implied WC=0 and DH ^ . HQ(-) contributes 
to a decrease in TEMS (an increase in the heat deficit) . 

The other properties of the pack remain unchanged if all 
other influencing factors remain negligible. 
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Figure 5: Simulation 
Negative 
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Snowmelt-No-Rain Positive HQ - Figure 6 is the algorithm for 
"Snowmelt-No-Rain, Positive Energy Flux", Thus, if TQ <1, 

WC >0, DH = and HQ(+) produces melt. The melt, MT, increases 

WC of the pack and produces excess melt^ EMT, if the WHC of 
the pack is exceeded. 

In the procees there is a reduction in DIPS due to the 
increment of melt and TEMS remains at QOc. If TQ >1, HC <0 

(DH^^O), and WC = 0; hence HQ'( + ) goes to satisfy DH and the 
remainder, if any, produces melt. 
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Figure 6i Simulation Algorithm for Snowmelt-No-Rain, Positive 
Energy Flux. 
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Snowmelt - with Rain, Negative HQ - Figure 7 shows the 
simulation scheme for "Snowmelt-with-Rain Negative Energy 
Flux. If TQ <l,it is implied that WC >0 and DH = . nQ(-) 
goes to freezing all or portion of the WC plus rain. L HF is 
released by the freezing liquid and changes TEMS accordingly 
If only a portion of the liquid water freezes, the remainder 
increases WC and drainage from the pack as EMT occurs if WHC 
is exceeded. 



If, on the other hand, TQ >1, it is implied that HC<0, WC=0 
and DH#0. HQ(-) freezes the rain and HC is reduced or 
increased depending on the amount of I^HF released. If only a 
portion of the rain freezes, HC - and WC is increased by 
the remainder of the rain. 
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Figure 7: Simulation Algorithm for Snowmelt- With- Rain, 
Negative Energy Flux. 
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SnQwmelt-With-Rain, Positive HQ - Figure 8 displays the 
simulation algorithm for "Snowmelt-With-Rain Positive Energy 
FIuk" 

For TQ <1, it is implied that WC >0, DH ^ and all of HQ(+] 
and HP produces melt. 

EMT results if AMT plus rain exceeds the WHC. There is a 
reduction in depth due to AMT and settlement. 

However, if TQ >1, the DH is satisfied before melt occurs. 
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Figure 8: Simulation Algorithm for Snowmelt-With-Rain, 
Positive Energy Flux, 
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Applicable Options i Data Input 

Table 1 outlines briefly the basic data requirement for the 
four different options formulated, Basic to all options are 
observed initial snowpack depth, density or water equivalent, 
precipitation, temperature and wind speed. In general, 

dewpoint temperature and vapO'Ur pressure are estimated from 
observed temperature and humidity with the use of peychro- 
metric tables (8) . 

Option 1 '. All of the available observed meteorological data 
are utxlized. In addition, computed estimates for potential 
evaporation according to a suitable empirical formula in- 
volving temperature and duration of sunshine or net radiation 

are ■used as supplementary data. 

Option 2 1 The requirements are similar to those of Option 1, 

except that the radiative energy supply is provided by 
observed net all-wave radiation. Initial TEMS, WC and TQ 
are assigned nominal estimates and WHC is estimated as a 
function of observed DENS. 

Option 3 : Estimates of shortwave radiation, (eq- {4))fare 

determined by using a nomograph of daily insolation for a 
given latitude, date and observed daily duration of sunshine 
(5)< The albedO' O'f snO'W pack is ass^umed to decrease exponen- 
tially accO'rding tO' the relationship: 

(24) ALB. At = ALB0.EXP(-At/KALB) 

where: 

ALB0 is the initial measurement or estimate of the 
pack's albedo; 

ALB is the computed albedo for time. At, subject to 

arbitrary limits of maximum (ALBMX) and minimium 
(ALBMN) values; and 

KALB is an empirical coefficient in time unit. 

For period with large amount of fresh snow, ALB is equated 
to ALBMX, Estimates of the longwave radiation (eq. (5)), are 
computed empirically in terms of air temperature (5) . 

Option 4 ; The basic available observed data of precipitation, 
air temperature and an index of wind speed are utilized. 
The remaining meteorological inputs are provided by computed 
estimates. Equation (4) provides for estimating HRSN where 
CHRS is determined from nomograph of direct insolation for a 
given latitude and date; and ALB is determined by eq.(M). 
Similarly, eg, (5) provides for estimating HRLN. In addition, 
the required values of CN are estimated from records of 
daily sunshine duration records. Values of the coefficients, 
KT, KG and KCL in eqs. (4) and (5) and estimates for snow 
cover depletion are guided by literature cited empirical 
relationships (5) . 



Table 1: Smranary of Data Requirements for Simulation Options 
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PR (rain) 
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TEMA (air temperature) 

HRS (shortwave radiation) 

HRL (longwave radiation) 

HRN (net all -wave 
radiation) 

ALB (albedo) 

HUM (relative humidity) 

VW (wind speed) 

SR/CN (sunshine duration) 

TEMD (dewpoint temperature) 

EA (vapour pressure) 

PET (potential evaporation) 

SOD (snow cover depletion) 
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Observed at site or extrapolated from nearby index station in the region, 
Psychrometric function (formulae/ table) . 
Empirical function (formulae/regression) , 
Regression estimate. 
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MODEL APPLI CAT I ON 

The basic input data requirements are the observed or 
estimated meteorological variables. These data provide for 
the sequential daily simulations. Snow survey data provide 
initial values for the snowpack condition and also aid in 
the overall watershed model calibration by way of comparing 
observed and simulated flows and snowpack variables. 

In addition to the data input options there are supple- 
mentary subroutines, formulated to perform computations for 
the interpolation of precipitation and other meteorological 
variables relative to snow course locations and to deter- 
mine, subsequently, weighted average inputs for given 
elevation zones. An inter station correlation coefficient 
and distance relation is used as a basis for the interpola- 
tion procedures . 

Option 2 is incorporated into the MOEHYDR model package 
which is tested in the Wilmot Creek IHD representative basin 
(see Part 2) . The following are brief discussions of the 
test results. 
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Test Results 

Figure 9 shows the simulated daily average snowpack depths. 

The solid trace represents the simulated values and the 
solid dots the maximum average value observed at a single snow 
course, while the open dots indicate observed basin average 
(12 snowcourses). values. Associated with each observed 
value is a dispersion of one standard deviation as indicated 
by the bar graph. 

In general, the simulations appear reasonable, in view of 

the large standard deviations associated with the observations. 
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Figure 9: Simulated vs Observed Daily Average Snowpack Depth, 
Wilmot Creek Basin: 1970-71, 
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Figure 10 shows the corresponding simulated snowpack water 
equivalent and density in comparison with observed values. 
The lower graph shows density and the upper shows water 
equivalent. 
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Figure 10: 



Observed vs. Simulated Basin Average Snowpack 
Water Equivalent and Density; Wilmot Creek Basin 
1970-71. 
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Figure 11 displays the corresponding simulated excess melt 
(lower graph) and evaporation/sublimation froin the snowpack 
and interception storage (upper graph) . 
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Figure 11: Simulated Daily Average Snowpack Excess Melt 

and Sublimation? Wilmot Creek Basin: 1970-71. 
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Figure 12 shows the comparison plots of the observed vs. 
simulated flows corresponding to the 1970-71 snow season in 
the Wilmot Creek basin. The melt rates, producing the 
excess melt to net basin input for flow synthesis, corres- 
pond satisfactorily with observed basin runoff responses. 

Some difficulties are experienced in matching simulated and 
observed flows during the periods of freeze-up and ice 
conditions. This is probably due to the irratic freeze-thaw 
cycles experienced in these parts of Ontario. Due to sudden 
decreases in diurnal air temperatures a portion or all of 
the excess melt drained from the pack is, usually, refrozen 
before it can be absorbed by the basin to produce runoff. 

These frozen ice layer conditions are subsequently remelted 
with the onset of favourable climatic conditions at a later 
date and contribute to total runoff. The existing structure 
of the hydrologic model do not simulate these phenomena. 
The daily excess melt as simulated, is immediately utilized 
by the model as net basin input to produce runoff. 
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Runoff; Wilmot Creek Basin i 1970-71, 
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SUMMARY AND CONCLUSION 

A snow accumulation and melt model is formulated, computer- 
ized and tested for its applicability in providing reliable 
simulation of excess melt for use in synthetizing and pre- 
dicting winter flows in selected Ontario drainage basins. 
Supplementary routines are formulated to accommodate the use 
of different simulation options depending on the availability 
of local climatic data. Based on calibration tests with the 
MQEHYDR-conceptual watershed model, the SNMELT submodel in 
application with data input Option 2, performed satisfac- 
torily in providing computation of index of net water input 
to the basin to produce streamflow. 

The model is formulated to relate to snow accumulation and 
melt experiences in Ontario. Nonetheless, it is sufficiently 
general to allow its application in other snowmelt regions. 
However, suitable local empirical relationships for average 
snowpack density, depth and selected climatic variables are 
required, and need to be carefully derived in order to 
ensure satisfactory application. The SNMELT model also 
lends itself to be used as a tool in simulation exercises 
for hydrologic design purposes. Because the initial con- 
ditions of the snowpack depth, water equivalent or density, 
can be easily updated from current observations, the model 
may be used as an auxiliary submodel in spring runoff 
forecast operations. 
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A, {AAP} 

A, (AMP) 
Ai, (A RAP) 

AIAARMP) 
AC A EI S 

AC A ESP 

AEIS 

A ESP 
AET 

ALB 

ALBMN 
A LB MX 
AMIS 
APR 
A VTEM 

B, (BAP) 

B, {BMP} 

Bi, B2, (BRA PI (BRA PI) 

fl/, B2, (BRMP), (BRMPI) 

BS 
CHRS 

CN 

CORLN 

CORLN0 

CW 

DDEP 



Notatifms 

^ regression conslanl in densily-deplh relalion lor accuniulation 

period, no rain 
= regression constant in density-depth relation for melt period, no 

rain 

= regression constant in density-depth relation for accumulation 

period, with rain 
= regression constant in density-deplh relation for melt period, with 

rain 
= aocumulation variable for actual sublimation/evaporation, irom 

interception storage 
= accumulation variable for actual sublimation/evaporation from 

snow pack 
= actual sublinialion/cvaporalion from interccplion storage 
= actual subliiiiation/cvaptvralion from soowpack 
= accumulation; variable for actual sublimalion/evapo'riilio'n fVoui 

interccplion sioragc and snowpack 
= albedo of snowpack 
■= assigned minimum vakjc for albedo 
= assigned maximum value for albedo 

= actual water-depth equivalent of melt in interceptioo storage 
= accum.ulated rain-on-snow for given period 
= variable for apparent average snowpack temperature 
= regression coeflkient in density-depth relation for accumulation 

period, no rain 

— regression coefficient in density-depth relation for melt period, no 
rain 

= regression coefficients in density-depth relation for accumulation 

period, with rain 
= regression coelllcienls in density-depth relation for melt period, 

with rain 

— assigned coeOicicnt for interception storage 
= clear sky insolation 

= cloud cover factor 

= intcr-stalion correlalion coelficient 

= assigned minimum value of intcr-statioo correlation coeOicient 

= specific heat of water 

= depth reduction in snowpack due to melt 
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DEN, (DENS) 
DEN I, DEN2 

DEN{NEW), {DEN0) 

DENW 

DH 

DIST 

D/Sr0 

'DISTEMP' 

'DfSTOMDTA' 

'DISrPPT 

or EM, {DTEMC), (DTEMF) 

DTEMIDZ 

EA 

ES 

EMA 

EMS 

EMT 

FAC 

FACT 

PARIS 

FWC 

HA 
MAC 

mc 

IfCOND 

IICONV 

IIGD 

HM 

HQ 

MR 

HP 

HRL 

HRLN 

HRN 

MRS 

MRSN 

HUM 

JS 

ISCM 

KALB 

KC 

KCL 

KCOND 

KCONV 

KF 

KGD 

KT 

KM. 



average density of snowpack 

average density of snowpack at the start and end of each time 
increment, respectively 
average density of fresh snow 
density of water 
increment of heat storage 
inter-station separation distance 
critical value of inter-station separation distance 
subroutine for temperature data input 
subroutine for other meteorological data input 
subroutine for precipitation data input 
change in average temperature of snowpack in "C, "F 
temperature gradient across ground- snow interface 
vapour pressure at given air temperature 
saturated vapour pressure at given air temperature 
emissivity of air 
emissivity of snow 
excess melt 

variable defining current quantity of energy, per degree tempera- 
ture of the pack's condition 

variable for defining temperature conversion from "F to °C 
fraction of drainage area contributing to interception storage 
potential value of liquid-water content that can be frozen by 
outgoing energy flux 

variable defining aerodynamic energy flux 
variable defining advective and conductive energy fluxes 
variable defining heat capacity of the pack 
heat of condensation/sublimation 
heat of convection 
heat of conduction 
heat equivalent of melt 

variable defining change in heat stored in the pack 
variable dcfming radiative energy Hux 
heat of advection of rain 
longwave radiation 
net longwave radiation 
variable defining net radiative energy flux 
shortwave radiation 
net shortwave radiation 
relative humidity 
increment of interception storage 
interception storage capacity 

empirical cocnicicnt for reduction in albedo with time 
empirical cloud ceiling factor for shortwave radiation 
empirical cloud ceiling factor for longwave radiation 
a transfer coefficient for heat of condensation 
a transfer coetficient for heat of convection 
forest cover density factor 

an assigned constant for ground heat conduction 
forest cover transmission coefficient 

empirical coefficient for change in average snowpack density due 
to compaction 
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LHF 

IJIFS 

LW 

MELV 

MGD 

MIS 

MIST 

Ml 

P, (PRECfP) 

PC AREA 

PEfS 

PESM 
PESMZ 

PESP, iPSE) 

PET 

PFWC 

PHE 

PN 

PR 

PI 

PRN 

PS 

PSN 

RMIS 

SR 

SIGMA 

SCD 

TEM, (TEMS) 

TEMl, TEM2 

TEMA 

TEMD 

TEM0 

TQ 

TQJ, TQ2 

TQ0 

VIS, (F/SO) 

yw 
wc 

WCh ^C2 

WC0 
VVEQ 
WEQh WEQ2 

WEQ0 

WGT 
WHC 
WHCI, WHC2 



WHC0 

At 

I 
J 
K 
H 



^ latent heat of fusion released by freezing rain 

= latent heat of fusion for ice 

= variable defining amount of liquid-water in the pack 

= mean elevation for elevation-zone 

= ground melt 

— variable defining incrcmenl of melt in interception storage 
= variable defining total melt in interception storage 

— variable defining increment of total melt 
= increment of total precipitation 

— elevation-zone area percentage 

= variable defining potential sublimation/evaporation from inter- 
ception storage 
= variable defining potential evaporation from soil moisture 

— variable defining potential evaporation from soil moisture for 
elevation-zone 

= variable defining potential sublimation/evaporation from snow- 
pack 

= potential evaporation 

= variable defining potential volume of liquid-water that can be 
frozen 

= variable defining potential heat of sublimation from snowpack 

= variable defining net rain remaining after freezing in the pack 

= increment of rain 

= variable defining net water input to the basin system 

= net increment of rain after interception storage 

= increment of snowfall 

^ net increment of snowfall after interception storage 

— variable defining remaining volume of melt in interception storage 
= hours of sunshine ratio 

= Stefan's constant 

= variable defining snow cover arcal depletion rate 

= variable defining average snowpack temperature 

= variables defining average temperature of the snowpack at the 

start and end of each time increment, respectively 
== average air temperature 
= average dewpoint temperature 
= variable defining average temperature of fresh snow 
= thermal quality of the snowpack 
= variables defining thermal quality of the snowpack for the start 

and end of each time increment, respectively 
= variable defining thermal quality offresh snow 
= variable defining current volume of precipitation in interception 

storage 
= wind speed 

= liquid-water content fraction of snowpack 
= variables defining liquid-water content fraction of snowpack at 

the start fifild end of each time increment, respectively 
= liquid-water content fraction of fresh snow 
= water equivalent of snowpack 
= variables defining water-equivalent of snowpack at the start and 

end of each time increment, respectively 
= water equivalent of fresh snow 

= distribution weight for snow course 

= liquid-water holding capacity of snowpack 

= variables defining liquid-water holding capacity at the start and 

end of each time increment, respectively 
= liquid-water holding capacity of fresh snow 

= time increment 

= I, 2, . . ., M = number of meteorological stations 

= 1,2,..., A^A^ = number of days of simulation 

= 1,2,.. ., MM = number of elevation-zones 

= 1, 2, . . ., m = number of snow courses per elevation-zone 
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COMPUTER DOCUMENTATION: SNMELT 

Table 2 gives a brief suiranary of some aspects of the com- 
puter documentation (11) that are directly associated with 
the use of the snowmelt model in relation to the MOEHYDR - 
hydrologic model package. Gomputer program listing for 
MODEL and SNMELT and description of the parameters/variables 
associated with SNMELT are also given in the Appendix, The 
remaining subroutines {I-NPUT, SUMMARY, etc.) and submodels 
(LTSM, DEW, STSM, etc.) are listed elsewhere (11). Alter- 
native user-oriented versions of these subroutines and 
submodel may be formulated and implemented with SNMELT 
instead. 

Sample input format for the model, basic data input and 
examples of the snowmelt outputs, {tabulatiQns and plots) 
are given in the Appendix, 
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TABLE 2: SUMMARY OF COMPUTER DOCUMENTATION FOR THE f^EHYDR-SNOWMELT MODEL 



NAME/USAGE 



MODEL (PAR, PREDY, NPAR, NOBS, 
ISTART, IGEN) 



PGM. CHARACTERISTICS 



SUBROUTINE 



PURPOSE 



INPUT, DEW, SNMELT LTSM, Read data files? call sub- 
STSM, SUMARY, ENDMON, QPTl, routines for modelling func^ 



LANGUAGE 

FORTRAN IV G 
Level 21 



0PT2, OPTIO 



tions and suinmaries 



MACHINE CORE USED 
/ 360-65 9.8K 



INPUT (ISTART) 



none 



Inputs parameter values, sets 
initial conditions and ac- 
cumulation variables 



- ditto - -ditto- 5.6K 



DEW (TEMA, RH, EA, TEMD) 



none 



Uses air temperature and 
relative humidity to calcu- 
late vapor pressure and dew- 
point temperature 



- ditto - -ditto- 5 . IK 



SNMELT (SNOW, ISTART) 



none 



Simulates the physical proces- 
ses of snow accumulation and 
melt 



- ditto - -ditto- 9 . IK 



LTSM (ISTART) 



none 



Simulates the watershed hydro- 
logic processes on a daily 
basis 



- ditto - -ditto- 6.9K 



STSM (ISTART) 



none 



Simulates the watershed hydro- 
logic processes on an hourly 
basis 



- ditto - -ditto- 11. 8K 



SUMARY (I, IIDAZ, ISTART, MONTH) DAYSTR 



Prints out summary of all 
significant variables on a 
daily basis 



- ditto - -ditto- 4.6K 



ENDMON (MONTH, ISTART) 



LINEP, 0PT3, 0PT4, 0PT5, 
0PT6, 0PT7, 0PT8, OPT 10 



Prints/plots end of month sum- F0RTRAN IV 
maries storages and flows G level 21 



-ditto- 



4.8E 
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TABLE 2 continued 



NAME/USAGE 



OPTl 



SUBROUTINE 



none 



PGM. CHARACTERISTICS 



PURPOSE 



LANGUAGE MACHINE CORE USED 



Tabulates monthly summary for F0RTRAN IV G /360-65 3.6K 



the different simulation 
blocks 



level 21 



OPT 2 



none 



Tabulates simulated daily 
flows 



- ditto - -ditto- 1.6K 



OPTIO (MONTH) 



0PT3(PRTMX, AMACC, MONTH) 



none 



none 



Prints out summary of model 
parameters 

Tabulates output of the snow- 
melt submodel 



- ditto - -ditto- 5.6K 



- ditto - -ditto- 3.3K 



0PT4 (PRTMX, AMACC, MONTH) 



none 



Tabulates daily summary of 
infiltration and soil moisture 
storages 



- ditto - -ditto- 3.6K 



OPTS (PRTMX, AMACC, MONTH) 



0PT6 (PRTMX, AMACC, MONTH) 



none 



none 



Tabulates monthly report of 
daily output from interflow 
and Baseflow storages 

Tabulates daily summary of the 
primary components of the 
simulation 



- ditto - -ditto- 3.8K 



- ditto - -ditto- 3.5K 



0PT7 (MONTH) 



none 



Plots of precipitation, soil 
moisture storage , evapotrans- 
piration, observed and simulated 
flows 



- ditto - -ditto- 2.5K 



QPT8 (MONTH) 



DAYSTR (DTRF, DRIFL, DRBF) 



LINEP (MONTH) 



none 



none 



none 



Plots output of SNMELT - depth, - ditto - 
water equivalent and density 

Metric conversion and temporary 
file storage 



Alternative plot for flows 
simulated vs. observed 



-ditto- 2 . IK 



2.7K 



- ditto - -ditto- 1.8K 
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NAME/DSAGE 



MAIN 



SUBROUTINE 



MODEL, OPTUM, LIMIT, 
WGT, CONV STATS, 
SURF AS 



PUKPOSE 



Controls the optimization and 
simulation procedures 



PGM. CHARACTERISTICS 



LANG. 



MACH- 



CORE 



7.7K 



OPTUM (lOPT, ISTART, 
SST, NOBS) 

LIMIT (TPAR, LPAR, NPAR, 
ISTART) 



WGT 



none 



Performs optimization 



Monitors the limits of the para- 
maters 



- ditto - - ditto - 0.9K 



- ditto - - ditto - 0.5K 



CONV (NPAR, TPAR, DR, B 
BOG) 



none 



Monitors the convergence of the 
parcimeters 



- ditto - - ditto - Q.7E 



STATS (NPAR, NOBS, SSE) 



none 



Calculates statistics of fit; 
observed vs simulated flows 



- ditto - - ditto - 1.5K 



SDRFAS (NPAR, NOBS, 
ICONV, lOPT) 



MODEL, OPTUM 



Calculates response surfaces for 
pairs of parameters 



- ditto - - ditto - 2.3K 
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PROGRAM LISTING 



' MODEL • 
' SNMELT ' 



V 



~rGT?7TriSi IV 5 LEVEL — rr 



0001 

0002 
0003 

oooa 

0006 
0007 

oooa 

0009 



SUfiROU 
CDMVOM 
CRMMOM 

rpTTusr 

COMrtON 
COMMON 

I, PAR IS 

COMMON 



COMMO 
COMMO 
COMMO 
COMMQ 




PAGE COOi 



0010 
0011 
0012 
0013 



0014 
0015 
00 1 6 

0017 



'OOIB 

0019 
0020 
Q02t 



0022 
0023 
002a 
0025 
00<i6 
0027 

0028 
0029 
0030 

iroTr 

0032 
0033 
OOjii 
003S 
0036 



GOMMO 

CnMHQ 

CO HMO 

liACUS 

CPHMO 




FfSHC 00 
00 
00 
000 
000 



JTZnTTTT 

oaseo 

04590 
04600 

oafeio 

0U620 



000 

000 

K,C 000 

fOACODOOO 

,OARIF00O 



04630 
04640 
04650 
04660 
04670 
04660 



v«r 



000 
OOO 
00 
000 

ooo 

000 



046«50 
4 7 00 
04710 
04720 
04730 
04740 



Tim 
000 
000 



04750 
04760 
04770 



CfJMHQM 
COf^MQN 
CGMMOM 
COHNQK) 
COMMQM 
OIMENS 



7CCC/ C0PH,CUSWI3W,CDPfc.T 
/CDaY/ ndav 

/CTPiNT/ rMPANT(10)|AHPANT(3) ,TmVIS 
/DMTRX/ PRTMJ(C33,31),AMNACCC2n 

/opm/ iourop(20) 

ION PARC20)rTPAR(20),OBSY(750)|PREDY(750) 



"00 57' 
0038 
0039 
0040 
004 1 
04? 

"0 04 3 
0044 
0045 
0046 
0047 
0048 



00004780 



DIMENS 
REAL I 
REAL K 
RFAL*a 
INTEGE 



ION HUURP{24J 

ION ■INDAT(750,2),DAT(750,9),OBSFRV(750),HOURPSC100,24) 

SC,ISCM,lS,K2,K3,Kt,MTiNN,SUNUP(l2),SUNDNC12} 

PA(?»PREUY,0BSY,8PAR,TPAR 
» STM0ATC400) 



00004790 

00004800 
4 810 
00004820 
00004630 



[H3~?"(rT 

9092 AMNACC 

NDAYsO 

" ■ 171=0 

DO 910 

910 PPEOY( 



■i JBRU = l,21 
CJB«0)sO,0 



J = r,MOftS 
J)=0,0 



00004850 
00004860 



4^ 



■>ORTRAK IV e LiVEL — FT 



'Mpiitt' 



ISTT^TfWf 



lE/25/a3 



PAGE 0002 



00ii9 
0050 
OObl 
0052 
~0'0S3~ 

005a 

0055 
0056 
0057 
0056 

iro"s^ 

060 



0061 
ITffbT 
0063 



006U 
"UObT" 
0066 
0067 



00 63 
0069 
0070 
00 71 
0072 



Trrs7 
oofiy 

0085 
0066 
0067 
00S3 



OOSi 
009 



irCISTART,NE,l JGO TO 9«0 
wf?lTEt6,20o2) 
2B02 F0KMAT£5 0FIRST PASS THROUGH THE MODEL') 
rtP I TE (6,6667) (IHON(KL) i KL= I * 12) 



onooafl70 

OOOOaSBO 
OO00aB90 



6667 
940 



fOi^HAT fix, 1215) ■ — --^ ■ ■ ■■ 

CONTI^JUE 
iMPfsS 
INPT1=5 
IhJpT2 = 5 

IKPTjsS 

I SIMS I ~" ^- -~ ■""""""" "" ' " ' '' _— ^ 

JCNTsl 

READ IN DATA TO BE INPUT If THE pPUT ROUTINE 

JCALI INPUT (ISTANT) 

lFtISTARr>E,l) GO TO 5 ~ ^''' ' ^^ "^~ 

READ tS, 10) IRH 

READ IN INFORMATION ON PIRlOD OF SIMyilTlON AND STOfiM DATES 

RE AO ( J NPT I 10) NSYO , NDAVS , LEAP 1 1 LEAP 2 



00004900 
0000^910 
0000^920 
OOOfiaQSO 

ooofia9ao 



OOOOU950 

ooooa96o 

OOO0U970 
0000a980 
0000«990 
00005000 

■flOlTOSOTO" 



' NEYOaWSYD+NUAYS-l 
10 FORMAT C20I a,/, 15(201 a,/)) 

RE AD CINPT # 1 ) N8T MO # CST MO A? CI ) , I « 1 # NSTmD ) 

READ IN OA|A ON SUNRISE AND SUNSET TIMES BY MONTH IN HOURS 



00005020 
00005030 
000050^0 
000050SO 



00005060 
00005070 
00005080 
00005090 
00005100 
00005110 



HEA0tINPT#20) (SUNUP (I)#Isl,12) 

READ£lNPTi20) tSUNON (n,I*i|12) 
20 FORMAT(12F5,0) 
5 CONTINUE 
ICNTsl 



DO PROCEDURE FOR NUNBER Or f^IVs iNPuT 



00005120 
00005130 

oooosiao 

00005150 
00005160 
00005260 




00005276 
000052S0 
00005290 



00005300 



00005310 
00005320 
00005330 
O0OOS3aO 
00005350 
00005360 

wuo"5r7Tr 

00005380 
00005390 
00005U00 
OOOOS^IO 
00005^20 



READ IN DATA RELl¥ANT TO OaV 

IFClRH,EQ,OrOR,ISTART,NE,l) GO TO 7661 

READ 1 1 NPT 1 ,331 ) ( INPAT ( IZl # JST) , JSTal ,2) , (DAT CIZI , JST) i JST^l #6) f 
I RH»DAT(IZl,9) 



■o7nro5^ir5T 

00005^40 
OOOQSaSG 



■< 



isr 



a*P? :-!-ft^^t" 



> 



f f-OHTRAfJ 14 G LEVEL JT 






,^- >. .^ 



si' *^iw: 



DATg 3 t?Q>7$ 



la/iS/ai 



PAGE 0001 



"^ 



'3 



91 

92 
0093 
OOVii 



1JT9T 
0096 
0097 
0099 
0099 

oion 
irrcT 

0102 
0103 

otoa 

0105 

OlOfa 
'TOT 
010« 
0109 
OHO 
01 11 
0112 
"OUT" 
OllU 
0115 
0116 
0117 
0118 



UTPr 

0120 
0121 
0122 
0123 

oia'i 



331 Fnf<M4T(ioX,l2, I2,2X,8F8g2J 

JMTM = rrvr)AT(TZl,l ) 
j:^Ay = !NDiT(I^l,2j 



DPET = 0AT(IZU3J 

CDPWsDPR 

CDPET=DP£T 
ADPR = ADPR + lJP*^ 



S 



TrrsvoviiiTjsisiuwf + usNOrt 

ADPETsAOPET+DPET 
SNOntPsDATClZt ,«) 
TEMA=0iTtlZlr5) 
HFa=DAT(IZl,fa) 
VWsOAT (IZ1,9) 



7&61 



CALL Dt>«CT£:MA,Whl,TtMD,|AJ 

0AT(IZU7)=TEMD 

0ATCIZl,8)=EA 

GO TU 76o2 

COMTINUf 

IFdSTART.NE, n GO TO 7669 



' ffOTrnTrpnT33 h inoa r c i ^ i . jit ) , jsts i , 23 » t D at 1 1 Z l , JST J , JST* l #9 J 

7669 JMTH=INDAT(lZlf U 

JDAY=INDAT(1Z1,2) 
"-*■' DPW = r)ATaZl,l) 

DSN0iN = DAT(IZU21 

0PET=DATtIZI,3) 



CDPWsO 

cnsNO'/- 

CDPtT= 
AnP(<rA 

A n s ^ w 

AOPETs 



0125 
0126 
0127 
0128 
0129 
0130 

0132 
0133 
013a 
0135 
0136 
"OUT 
0138 
0139 

oiao 
oiai 
1 a 2 



1 a a 

.0145 

oia7 

D14B 



£T 

R+DPR 

DSnOw+DSNOW 

PET+DPtT 



TtMAsD 
WFucHA 

EAsDAT 
Vw=DAT 



"TTO" 
7662 



7665 



981 

332 

30 



CCNTIhJ 
IPCIST 
REARCI 
OflS»=0 
IF(TST 

CCNTIM 
ICM = I 
FORMAT 
FOKMAT 

DO ao 



^nTTTzrj 

(lZt/5) 

IZt#6) 

{IZ1,7) 

21*6] 

Zl,9) 

^"X,12,I2,2;i,7F8,2,2F4,0j 



40 
50 



"TTTTDTT 
c n N T I N 

MONTHS 
SUNRIS 
SUNSET 
IFCIOA 



T.nE.U go to 7665 
T2,332) OBSERVdZl) 

E«v(izn 

r.NE.n GO TO 981 

rrn"=u¥^i? ' 

T + l 

F10,2) 
Flo. 21 
1 «12 

"LErrnmrrrrTTTjirTirTr 



UMjP(MONTH) 
UNONCHIQNTH) 
EQ.STMQATCISTM)) GO TO 60 



00005500 



00005530 



0~G 0554 0" 

00005550 

00005560 

00005570 

OO00S56O 

00005590 

Tn)"o"D"s^y<r 

000n561f> 
00005620 
00005630 
000056*40 
00005650 






s.Jbi 



FORIRAM IV G uEvEL IT 



rtOOEL 



OAT£ = 77fi75 



U/2S/U3' 



PA5C oooa 



^ 



01U9 

01 so 

Olbl 
0152 

OlSa 
0tS5 
015b 
01S7 
01i8 



IF(TSTAST,eQ,!) GO TO 230 

PET = tJPE.T*AkEA*DELTAT 
NSTEPsl/DtLTAT 



OOOObfcbft 
0(^005670 
00006680 
0O00b6<>0 



IROJTE=0 

SKOi«snSNO"*AREA 

PR30PW*AHfcA 

on Po Jsi.^STfep 

IF(SNOW,Lt;,O,O,AN[>,D£Pl,LE,0,O) GO TO 75 
C^LL SNMtLTCS^Oi^/ISTABTl 



1)0005700 
00005710 
00005720 
00005730 
00005740 
00005750 



IPuUTtsl "■ ~ 

75 COriTlNUE 

CALL LTSM 
80 CONTISUE 

CALL SUMARYCliIDAZ,ISTART,M0NTH3 

GO TO 130 



OlbO 
0161 
0162 
0163 
OlbU 



"OTbT 
0166 
0167 
OlbS 

0170 
DT7T" 
0172 
0173 
017a 

0175 
0176 



■ST 



766fa 
7fa67 



DFLTAr=flinTEi; ■ • ■ 

PETsDPtT*ARE&/( CSUNDN tMON|M )-SUNUP (HONTH J ) /DELT AT ) 

MSTEPs 1 /DELTAT 

IF(ISTART,NE,1 J GO TO 7666 

KEAD(If>JPT3i I to) (HaURPS(ISTM,IJK),IJK=l,2aj 

on 7bb7 1JK=1,2« 



O0UO57bO 
00005770 
00005780 
00005790 

000056 10 
0OOO"5g2T 
00005830 
000058^0 



H0lJ9P(lJi<)3MuU«PSfISTH,IJK) 

110 FORMAT (16X,6F8,2,/,lfeXt8F8, 2, /,16X,8F8,2) 
ISTMsISTM+l 

IFCISTART.EQ.n GO TO 230 
00 100 IH0URsl,2a 
PRsHOURP(IhOU»)*AREA 



irrrr 

0178 
0179 
0180 
0161 
0162 



irTST 
018a 

01 as 

0166 
0187 

oiee 



Do 120 Jsl,NSTEP 

CALL STSM(ISTART) 
120 CONTINUE 
100 CONTIMUE 

CALL SUMARYC2,IDAZ/ISTART, MONTH) 
130 DO lao 1=1 #12 



IF(IDA2,tU.lM0JytnJ GO TO 150 

100 CONTINUE 

GO TO 230 
150 CALL END«D^JCMONTH,ISTART) 
DO 9091 JBRDal ,21 
9091 AVNACCtJ0RD)aO,O 



00005670 
00005880 
00005890 
00005900 
00 059 to 



00005920 
00005930 
00005940 
00005950 

00005970 

00005990 
00006000 
00006010 



0189 
0190 
0191 
0192 
0193 

0i9a 



ND A Y = ' 

299 CONTINUE 
230 CONTINUE 

DO 7070 1X1 = 1, ni 
7070 PRtDY(IXU=PREOY(IXl J*(AREA/2a.) 

IF(LEAP2,EQ,0) GO TO 8621 



ITYfT 
0196 
0197 
0198 
0199 
0?00 



"DO B622 IZVa2,l2 ' 

8622 rMON(IZY)siMON(IZY)-l 

8621 CONTINUE 

rF((ISTART.EQ,2),AN0,(IDUTOPCl).E0,in CALL OPTl 
IF((ISTART,E«,2),ANO,{IOUToP(23.EO.n) CALL 0PT2 
IF((ISTART,EQ,2),ANO,(IOUTOP(lO),EQ.l))CAiL OPT 10 
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END 
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COhhOW y'CI^iCS/ VIS, TP, TPn , iSCh, T ATrSTSC^TTITtS , aCpF tS, ACDTIS , iSC OOOOOlfTo" 
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COMMON /C 
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PhT,TE^47" rEMD,weur,WPT7t5ENi,RP4,K,C < 
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RIOn fOEGREES F) 
EES F) 

SNOWPACK (ACRE-INCHES) 

XGfc (ACRE-nrChES) 

STORAGE CACPE-INCHES) 
ITY CACRE-It^Cl-eSJ 
I.S, (ACRE-INCHES) 
(ACRE-INCHES/DAY) 
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00009500 
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00009520 

00009530 
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002H 
0029 
0030 
0031 

"D733r 
0033 
003« 
0035 
003b 
0037 
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0039 
0040 

ooai 
ooa? 
ooaj 



"DTTinr 

00«5 
0OU& 
00C7 

00«9 

"oiTSir 

0051 
0052 
0055 
005U 
055 



0057 
005B 
0059 
0060 

0061 



DVIS 
AfAET 
ACt)VIS 
ACIS 

TTtTTT" 



-I*JCRE"£NT OF VOLUME TO I.S C ACRE-INCHES 3 
-ACCU'^ULATED ACTUAL E.T, FROH ALL SnURCES UCRE- 
•ACCUMULATEO INCREHINTS OF vOLUf^E TO I.S. (ACRE* 
-ACCUMULATED VOLUME_ TO^ I,$, t ACRiE-InCHES ) 
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INCHf S) 
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00009610 
00009620 
00009630 






AT 



OF 



-^OTAT 
OFPl "DtPTH OF SNOwPaCk 
ISO=ISCM«A«EA*OELTAT*FAHIS 
hTsPFU 
TPMlsTEMO 

IRAIO ^^ 

IF"("p"ft7^t7r7T~"i77r:ri 

TPf*ECtPR+3N0'-N 

MCO^O = Hl*(EA«6.n)*V.»*203,? 

hCOMV=B2*(rEHA-3^,)*VW*203,2 

wHCUO,03 + 0,03*OtM 

TQlsl.0 -CTEMi,3g,0j*(S. Q/9,0)*(l . 0/160,01 



FROM I.S, (ACRE-INCHES) 
PERIOD (INCH If) 
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00009680 
00009690 
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00009700 
00009710 
0009720 
00009730 
000097^0 
00009750 
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IFCTPRtC.LE.O.O) GO TO 
IFCSNOw.LE.O.) GO TO 5 

™ IRAaIRA*2 
GO TO 39 

HOI AEISsO, 



aOl 



FTTsTT; 

SNO^NeO I 

GO 70 a) la 
39 CCNTI^JUE 

AEIS=0, 

OENOaO.lQ 

wHCO=0,03 

T^O=1,0 

TEHn=TEMA 

DeP0=SNOH*0ELTAT/CDENO*AREAl 

PKN=TPkEC 



O0009T6T5 

00009770 

00009780 

00009790 

00009800 

00009A10 

000098?15~" 

00009830 

00009B«0 

00009850 

00009860 

00009870 

■i5"(ro"0"9"afl6 

00009690 
00009900 
00009910 
00009Q20 
00009930 



n 



TS^By^^H*OtLTAT 
VIS=VISD+IS 
GO TO 27 
P^n'sTPREC 
IS=95*PWN*DELTAT 

vis=visn*is 



0062 
0063 
006^ 
0065 
0066 
0067 

"D'o fair 

0069 
0070 
0071 
0072 
0073 



/u 
0075 
0076 
0077 
0078 
0079 



UFPiJsxrro" 

WEiJsO.O 
DENn=0,0 
wCO=0,0 
TfJOsO. 

" sTlNuE 



-VISO 



27 COKTINUE 

IF(VIS,GT.ISO) VIS=IS0 

PEISsPET*FaRIS 
"" ' ' IFCIRA,EG,n GO TO 844 

IF(IRA,EQ,3) GO TO 802 
SKjOI«N = SNUh»(IS/{0ELTAT) ) 
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00009950 
00009960 
00009970 
00009980 
00009990 
TOTrOODT 
OOOIOOIO 
00010020 
00010030 

oooiooao 

00010050 
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0001 
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OOyl 
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T)060 
0070 
0080 
C090 
0100 

Olio 



p^sTTrrr 

60 TO 803 
ia« PK' = PR«(IS/DELTAT) 

SNOWKJsO, 
GO TO 603 
S02 PN=PR-(IS/C0ELTAT*2,)) 
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0001 
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S^Ow^' = S^^U'N-C^S/(DELTAT*^,) ) 
803 CONTINUE 

If (TE^iD,,GT,32,) go to 101 
MISTsn, 



00010180 
OOOlOl^iO 
001020 
0010210 
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Pht=rCOKD 


00010220 N 
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PSEsAflSCPHE)/203,2 
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oooioaao .. 
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IFCPEIS.GT.VIS) GO TO 102 


00010250 
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00010260 
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iDTAEIS = DTAEIS*AElS 
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V 


0137 




. DtN=OENl*Cl .♦K*DENO*DEPO) 
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T)T5T" 
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ii\iU CONTINUE 

IFCISTAR 

WRITE(6, 

5b7tt COMTIWUE 

311 F0»M4T(l 



T.NE.aj GO TO 5b7R 

3t 1 )PP,SNn-(,DEPlinEMl*TEMl,PN,S'^DWN#&ElS,HIST,VlS 



t TEMl = 
' MIST 



H0.< OUTPUT FOK INTEHCEPIION STORAGE FOR SNO<^MfcLT ROUTINE 
-' ',^7,0,' SWUW s ',F7,0, ' DIP I ■ '",F6,2,» OEWl s »,Pb,« 
•,Fb,2,/,» PN = ',F7,0,» S»JOwN s i,F7,0,' AEIS = ',F7,0, 

a 'iF7,0,t VIS = *,F7,05 
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AESP -ACTUAL EVAPOTRANSPIRATldN FROM THE SNOWPACK (ACRE-INCHE 
ACAESP -ACCUMULATED ACTUAL E V APOTR ANSPIRaT ION FlROM SNOwPACK (AC 
ACAET -ACCUMULATED ACTUAL E,T, FROM ALL SOURCES f ACRE-rNCHES) 
PFSP -POTENTIAL E,T, FROM SNOWPACK t ACRE-INCHES ) 
^mn? iFRACTlUM" 



*000l0''l6" 

00010930 

S00010930 

-oooioQao 

00010950 
00010960 



"C mn? iFRACTlUM OF 

E EMT -EXCESS MELT 
Z ACEHT -ACCUMULATED 

4 CONTINUE 

IF (TEH0,GT,32,) GO TO 
PHEshCONtD 
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113 



000109715" 
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00010990 

OOOllOOO 

OOOUOIO 

0001 1020 



€t 
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lia 



T-SEsiBS(PHEJ>203,2 
PESP=PSE*(1,-FARIS) 
IFCwEQl .GT.O.OOOn GO TO U^J 
WEU1=0, 
GO TO 115 
AFSP=(DEN1**HF)*PFSP 
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TmXfTPtAktAJ.GE, (PET-AtlSJl A^SPs fPf T»A|IS)/AREA 

IFCAESP.GT.wtOl) GO TO 115 
AC AESP = AC AESP + AESP* ARE A 
AC AET= ACAET ♦AESP* ARE A 
wEUls-^EQl-AESP 
DEP1=0EPI-AESP/0EN1 



Pl=PEPl-AESP/OENl 

x^JHEs(PHE*{l,-FAHIS))/^0i.2 
0TEM = )(PHE/(2,5a*rtEQl*C)*(9,/5,3 
TEMUTEMl+DTEM 
IF(T£M1,GT,0,] GO TO 48 

Tf^l=0, 
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Tf^l=0 
HG'JsO^ 

GO TO 121 

AESPsi^EQl 
k^EQUO.O 
DE?1=0,0 
DFnUO.O 
— -i" 1 - A — n 



TC 1 = , 
TOl=0, 
WhCl=0 ,0 

acaesp=acaesp+aesp *area 
acaet=acaet*aesp *area 

EMT=0,0 
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DTP5PR*4EIS+AFSP*AREA 
GO TO 991 
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0200 
0201 
0202 



020« 
0205 
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JOB time/CPU time 

FORTRMI I¥ G LEVEL 21 
MACHINE-/ 3 60 -65 
O.Sp - ASP 

The total transaction (job) time of the MOEHYDR-Model 
package is relatively small and very reasonable in relation 
to the size of the program package (core used - 220k) . The 
Job time for a straight simulation run with no optiraization 
is significantly less than an application with the inclusion 
of automatic optimization. Table 1 shows examples of 
transaction times. 



Table 1 Examples of CPU/Job Times 



JOB DESCRJPTIQN CPU TIME JOB TIME 

5 months daily simulation of one rvm. 4*2 sec. 20.1 sec 

no optimization with SNMELT basic 

output . 

Options s 0PT3, 0PT7, 0PT8 plus STATS 

7 roonths claily simulation of 25 optim- 1 min. 44 sec, 4 min. 56 sec, 
ization iterations with all Options 
plus STATS 



: ^ - 



INPUT/OUTPUT SPECIFICATIONS 
Model Control Cards 

Card 

#1. c.c. 1-4 - IRH (type of data set) 

= data set contains dewpoint temperatures, 
TEMD, and vapour pressures, EA. 

= 1 data set contains relative humidity 
instead of TEMD and EA 

#2. c.c. 1-4 - NSYD {starting day - Julian - of simulation) 
5-8 - NDAYS (number of days of simulation) 
9-12 - LEAPl {= year 1 not a leap year, = 1 year 1 
is a leap year) 
13-16 - LEAP2 (= year 2 not a leap year; = 1 year 2 
is a leap year) 

#3. c.c. 1-4 - NSTM (number of storm days to be simulated) 
5-8 - STMDAT(l) (day of year of 1st. storm date) 
etc. etc. 
69-72 - STMDAJ(18) (day of year of 18th storm date) 
etc. 

#4. c.c. 1-5 - SUNUP (1) (sunrise time in hours for Jan.) 
etc. etc. 

c.c. 56-60 - SUNUP (12) (sunrise time in hours for Dec.) 

#5. c.c. 1-5 - SUNDN(l) (s\mset time in hours for Jan.) 
etc. etc. 
c.c. 56-60 - SUNDN(12) (sunset time in hours for Dec.) 
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Data Input Cards 

If IRH = 0, Card #1. format is used 
If IRH = 1, Card #la, format is used 



- YR (year) 

- JMTH (montlil 

- JDAY (day) 

- DPR (daily rainfall .. in) 

- DSNOW (daily snowfall .. in) 

- DPET (daily PET estimates . . in) 

- SNODEP (daily snowpack depth .. in) (- Optional) 

- TEMA (daily average air temperature .. °F) 

57-64 - RP4 (daily net all-wave radiation -ly) 

65-72 - TEMD (daily dewpoint temperature ..'^F) 

73-76 - EA (daily vapour pressure ...rab) 

77-80 - VW (daily average wind speed . . mph) 



Card 




#1. 


C.C. 8-10 




11-12 




13-14 




17-24 




25-32 




33-40 




41-48 




49-56 



#l(a) c.c, 



8-10 


- 


YR 


11-12 


- 


JMTH 


13-14 


- 


JDAY 


17-24 


- 


DPR 


25-32 


- 


DSNOW 


33-40 


- 


DPET 


41-48 


- 


SNODEP 


49-56 


- 


TEMA 


57-64 


- 


RF4 


65-72 


- 


RH (daily average relative humidity 


73-80 


- 


VW 



%) 



#2. c.c. 1-10 - OBSR (observed daily flow cfs) 
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Optimization Control Cards 

Card 

#1. c.c. 1-5 - No. of data set 

#2. c.c. 1-80 - Title 

#3. c.c. 1-5 - No, of parameter tO' be optimized 
6-10' - No. of observations 
11-15 - Maximuni No_ of iterations 
16-20 - Objective funotion 

= 1 (Z squared residuals) 

= 2 (weighted I squared residuals) 

=3 (E predictions) 

- test level o^f convergence O'f parameter 

- standard value for decrementors 

- = 1 - optimize STSM par^ameters 
= 2 - optimize LTSM parameters 





c.c. 21-25 




26-301 




31-35 


#4. 


c.c, 1-10, 


#5. 


c.c. 1-5 



,1-20 ... /61-70 - initial estimate of parameter 1/2/ /7 

- = 0' actual data set used 
= 1 generated data set used 
6-10 - ■- 0/=l/=2/or=:3 - analysis with no response surface/with 

response surface/only response surface/or only prediction 

#6. c.c. 1-10/11-20/.. ../61- 70 - initial increment to perturb 

parameter 1/2/.../7 
#7. c.c, 1-10/11-20/. . ./61-70 - positive decrement constant for 

parameter 1/2/. , ./7 
#8'. - indicates output option (- not required; = 1 required) 

c.c. 1/3/5/7/9/11/13/15/17/19/21/23/39/ - OPTl/0'PT2/OiFT3/OPT4/OPT5/ 
OPT6/OPT7/OPT8/LINlP/OPTlO/SU'MMflfiY/ENDMON/HOU:RIiY 

#9. c.c, 1-2 - basin name 
21-40 - gauge name 



Model Parameters Cards 



Card 

#1. c.c. 1 - Code = "A" 

2-11 - DAYDEL (daily - time unit) 

12-21 - HRDEL (hourly-time unit) 

22-31 - AREA (basin area) 

#2. c.c. 1 - Code = "B" 

2-11 - WEQl (initial snowpack water equivalent ... in) 

12-21 - DEPl ■• •" depth ... in) 

22-31 - DENl '" * density ) 

32-41 - WCl " « water content ) 

42-51 - TQl """ '•* thermal quality ) 

52-61 - IC (aesigned coefficient - pack compactiori/settleiient) 

62-71 - C (specific heat of snow) 
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■ t 13 c.c. 1 - Code = "B" 

2-11 - BS (assigned coefficient for interception storage) 

*■ 12-21 - Bl (coefficient - heat of condensation) 

22-31 - B2 ( " - heat of convection) 

■ 32-41 - MISTl (initial melt in interception storage, IS) 

xy; , 42-51 - BF (coefficient for snowpack evaporation) 

1^^ ^ |4 c,c. 1 - Code ^ "B" 

t*'"; 2-11 - AM (density-depth regress, const. -Ace. period) 

^ *" 12-21 - BAP ( " " coeff. -Ace. " ) 

V ; ' 22-31 - AMP ( " " const. -Melt " ) 

^'" 32-41 - BMP ( " " coeff. -Melt " ) 

'^.j- 42-51 - ARAP( " " const. -Acc/rain " ) 

p; ■ 52-61 - BRAP( " •• coef. -Acc/rain " ) 

62-71 - BRAP1( " " coeff .-Acc/rain " ) 

fS c.c. 1 - Code = "B" 

f » 2-11 - ARMP (density-depth regress. const. -melt/rain period) 

|?l 12-21 - BRMP ( " " regress. coeff. -melt/rain " ) 

P^ 22-31 - BRMP1( " " regress. coef. -melt/rain " ) 

^' 

^- 

J: #6 c.c, 1 - Code ^ "C 



H/^ii 



2-11 - ISCM (interception storage (IS) capacity 

12-21 - PARIS (drainage area fraction of IS) 

22-31 - VIS (initial storage in IS ... in) 

if e.e. 1 - Code = "D" 

2-11 - DSCM (depression storage (DS) capacity) 

12-21 - FARDS (drainage area fraction of DS) 

22-31 - DVS (initial storage in DS) 

m c.c. 1 - Code ^ "E" 

2-3 - NLAYER (no. of layers of soil profile) 

5-5 - lYEARB (beginning year of simulation) 

8-9 - lYEARE (final year of simulation) 

ft c.c. 1 - Code = "E" 

2-11 - CP (Coeff. of infiltration - LTSM) 

12-21 - EP (Coeff. of infiltration - LTSM) 

22-31 - EM (Coeff. of infiltration - LTSM) 

#10 c.c. 1 - Code - "E" 

2-11 - FCAP(l) (field capacity of soil layer 1) 

12-21 - FCAP(2) (field capacity of soil layer 2) 

22-31 - FCAP(3) (field capacity of soil layer 3) 
etc. etc. 

#12 c.c. 1 - Code - "E" 

2-11 - ANTSM(l) (antecedent soil moisture - layer 1) 

12-21 - ANTSM(2) (antecedent soil moisture - layer 2) 

22-31 - ANTSM(3) (antecedent soil moisture - layer 3) 
etc. etc. 



/»,i! 
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#13 Q.Q. 1 - Code = "F" 

l; - 2-11 - FCAPIN(FCAP of interflow zone (IZ) 

l^f 12-21 - CI (critical volume ratio for IZ) 

^^ 22-31 - K2 (routing lag-const, for IZ) 

32-41 - TST (total storage in IZ) 

^ 42-51 - SHC (percolation rate for IZ) 

52-61 - RIFL0 (initial routed interflow) 

.) ■ 

#14 G.c. 1 - Code = "G" 

2-11 - K3 (routing lag-const for baseflow) 

\ 12-21 - RBFI^ (initial routed baseflow) 

f 22-31 - BFSM (storage upper limit of baseflow) 

32-41 - BP0 (initial baseflow storage) 



#15 c.c. 1 - Code ^ "H" 

2-11 - D0 (initial depth of over land flow) 

12-21 - OLF0 (initial overland flow) 

22-31 - AL (overland flow const.) 

32-41 - NN (overland flow coeff.) 



#16 c.c. 1 - Code = "H" 

2-11 - RSE0 (initial routed surface runoff) 

12-21 - USR0 (initial unrouted surface runoff) 

22-31 - VSR0 (initial vol. routed surface runoff) 



• #17 c.c. 1 ^ Code = "I"^ 

\ 2-11 - B (evapotranspiration coeff. - soil layers) 

12-21 - FINT (steady-state infiltration rate) 

22-31 - WP(1) (Wilting point-soil layer 1) 

I 32-41 - WP (2) (Wilting point- soil layer 2) 

42-51 - WP(3) (Wilting point -soil layer 3) 
etc,, etc. 

[ #18 c.c. 1 - Code = "J" 

t:- # 2-3 - M (No. of time-area-histogram (TAH) elements) 

U 

t #19 c,c, 1 - Code = "J" 

; ^ 2-11 - Kl (routing lag- const- for channel flow) 

ii ^ 12-21 - UOIiF0 (initial unrouted overland flow) 

t' ' 22-31 - TARHJS (1) (element 1 of TAH) 

[ 32-41 - TARHIS (2) (element 2 of TAH) 

I etc. etc. 

'.. 62-71 - TARHIS (5) (element 5 of TAH) 

\ 
I 

'■ #20 c.c. 2-11 - TARHIS (6) (element 6 of TAH) 

[ etCp etc. 

f ^ TARHIS (M) (element M of TAH) 

i #21 c,C, 1 - Code = "K" 

2-11 - TRANS (index of aquifer transmissivity) 

: 12-21 - ALITL (index of baseflow horizontal flow path) 

22-31 - ALARGL (index of baseflow longitudinal flow path) 

. 32-41 - HMIN (index of min. avg. depth of water table) 

J, 42-51 - HACT0 (initial depth of index of avg. watertable depth) 

52-61 - CS (index of avg. aquifer coeff. of storage) 
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- Part 2 « 

MOEHYDR-Snowmelt Model Application in 
Ontario Drainage Basins 

L. A, Logan 



JIBSTRACT 



The mathematical snow accumulation and melt model 
developed in the Water Modelling Section of the Ontario 
Ministry of the Environment is implemented as a comple- 
mentary submodel of a conceptual watershed hydrologic model 
for streamflow prediGtion and simulation under winter 
conditions in Ontario, The reproducibility of the model 
appears to be satisfactory as demonstrated by its ability to 
adequately simulate snow accumulation and melt experiences 
within a given IHD representative basin under varying and 
different winter conditions. The application is extended to 
other selected drainage basins located in different physio- 
graphic regions within the larger Lake Ontario drainage 
basin- The results are equally satisfactory and encourag- 
ing. The predictive ability of the model for possible use 
as an aid to flood forecasting models is demonstrated by 
simulating possible maximum spring runoffs from maximum 
snowpack water equivalent storage under exposure to critical 
sequences of climatic conditions favourable to snowntielt. 
The sensitivity of the snowmelt model outputs (accumulation 
and ablation), to errors in precipitation inputs, is shown 
to be highly significant for both underestimation and over- 
estimation of these inputs. 
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MOEHYpR-Snowmelt Model Application iii 
Ontario Drainage Basins 



L. A. Logan 

INTRODUCTION 

Several different approaches to the modelling and 
estimation of snowmelt for the synthesis of winter stream- 
flow have been attempted ( 1 , 2 , 3 , 5 , 6 , 10 , 13 ) ^ . The degree-day 
method has been the most widely used (5, 13) in computing 
snowmelts, because of its simplicity and ease of application 
in relation to data availability. Other approaches utilize 
generalized empirical snowmelt equations, which consider 
temperature as the primary heat index of melt in addition to 
some aspects of the physical properties of the snowpack 
(3,4,5), The more advanced approaches consider the theories 
of heat exchange by energy balance and aerodynamic princi- 
ples and the thermal properties of the snowpack in the 
modelling of the snowmelting processes (1,2,5,6,10), 

The purpose of this paper is to demonstrate the appli- 
cation of the MOEHYDR- snowmelt model under real time test 
situations in Ontario. The test applications are related to 
reproducing historical snowmelt events and snowmelt runoff 
predictions under superimposed critical sequences of meteor- 
ological variables. Also, the output sensitivity in the 
simulation of accumulation and melt to sampling errors in 
the estimates of daily precipitation inputs is also 
demonstrated . 
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SNOWMELT MODEL 

Snow accumulation, evaporation or sublimation and the snow- 
melting processes are considered as separate but linked 
events in the overall computational scheme. 

Figure 1 shows the general schematic flow chart relating to 
snow accumulation and ablation. 

Indicated are precipitation input as rain and/or snow, 
provision for interception storage, evaluation of the energy 
fluxes, accumulation and compaction of the pack and provisions 
for evaporation or sublimation from the interception storage 
and the snowpack. The scheme continues with computation of 
the potential energy equivalent of melt, snowmelt with 
negative or positive energy flux with or without rain on 
snow. Excess melt is subsequently drained from the pack to 
provide net water input to the basin for the synthesis of 
streamf low. 
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Figure 1: Simulation Scheme for Snowmelt Model. 






MODEL APPLICATION 

The snowmelt model is tested in the Wilmot Creek and 
Oakville Creek IHD Representative Drainage Basins for its 
applicability and suitability as a complement to the 
MOEHYDR-watershed hydrologic system model . Other drainage 
basins located within the larger Lake Ontario drainage basin 
were also selected as test basins. The physical energy 
equations employed are for snowmelt at a point; however, 
provisions are made for the use of average climatic data 
input of several point measurements in an attempt to account 
for spatial variability of these input data. Simulations on 
an areal elevation zone basis may also attempt to represent 
areal accumulation and melt rates. 

Measured and/or estimated meteorological data provide 
the sequential input requirement. Measured snow course data 
provide for the initialization of the snowpack condition for 
simulation and aid in the hydrologic system model calibration 
through the comparison of observed vs simulated flows and 
observed vs simulated snowpack conditions. 

The results to be discussed are related primarily to 
the model application with data input option 2(10). 

Dr a ina ge Ba s ins 

Figure 2 shows the general locations of the drainage 
basins used in the test applications. 

Table 1 summarizes the primary physical features of 
each drainage basin. The physiography of the Wilmot Creek 
basin comprises mainly interlobate kame moraine and till 
plain; the Oakville Creek basin is mainly till and clay 
plains and borders on the escarpment, while the Twenty Mile 
Creek and Welland River basins are mainly clay plain. The 
vegetation of these basins comprises croplands and pastures 
with sparsely distributed forested areas in the higher 
elevations , 

Input Data 

The basic input data for the simulation model are: (a) 
initial condition of the snowpack, (b) daily sequences of 
meteorological variables for the period of simulation, (b) 
(c) topographical and landuse data for assigning elevation 
zones and for determining areal distribution of forest cover 
and open areas, respectively and (d) snow cover depletion 
data for proportioning runoff from the snowmelt contributing 
areas . 



Direct observations of the snowpack provide for initial 
DEPS, DENS, TEMS, WC, TQ and WHC, Usually, measured values 
of the TQ and WHC are not available, initial values may be 
computed as a function of initial TEMS and DENS, respectively, 
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Figure 2 : - Location of Drainage Basins 
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Table 1: Physical Features of Selected Drainage Basins in Ontario 







Location 
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Figure 3 shows, as examples, variations in mean daily net all^wave 
radiation, air temper'ature and dewpoint temperature observed 
during the 1970-71 snow season in the Wilmot Creek basin. 
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Figure 3: Observed Average Daily Net All-wave Radiation, Air 
Temperature and Dewpoint Temperatures in the Wilmot 
Creek Basin: 1970-71 Snow Season. 
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Figure 4 shO'WS the corresponding daily variations in computed 
mean daily vapour pressure, wind speed and estimated potential 
evaporation for the same period (1970-71) . 
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Figure 4:' Average Daily Vapour Pressure Wind rSpeed and 
Estimated Potential Evaporation; Wilmot Creek 
Basin: 1970-71 Snow Season 



Historical Simulation 

Table 2 shows a statistical summary of the simulation 
for historical streamflows during the months of December to 
April, for four winters in the Wilmot Creek basin. The 
results are abstracts from the simulation runs for an over- 
all calibration of the watershed model. Comparisons between 
observed and simulated mean flows appear reasonably good for' 
the majority of the months. There are, however, some 
significant differences. Some of these differences could be 
explained by errors introduced into the interpretation of 
recorded flows by doubtful corrections for freeze-up and ice 
conditions for extended periods during several months. 
Other errors introduced are due possibly to unsatisfactory 
calibration of the watershed model. Also, the irratic 
short time freeie-thaw cycles that are frequently experienced 
in this region of Ontario cannot be adequately represented 
by the model. The results of the hydrologic model calibra- 
tion tests, however, indicate general acceptability, as the 
mean annual flows and the standard deviations of the monthly 
flows, for the corresponding observed and simulated flows, 
in most cases, are not statistically significantly different. 

Figures 5,1, 5.2, 5,3 and 5.4 show other results of the 
simulated snowpack depth, density or water equivalent and 
excess melt for the winter periods of 19 69-70, 1971-72 and 
1972=73. On the whole, the simulated snowpack conditions 
show good comparison with the observed values, in view of 
the large sampling errors that are associated with the 
spatial variation and distribution of the snowcover in the 
basin. 

Figures 6,1, 6.2, 7.1, and 7.2 show additional examples of 
the results of the model application in the other test 
basins. The simulations for the Oakville Creek basin, in 
general, appear to be satisfactory. The averages of eight 
snowcourses appear to be poor indices for the basin snow- 
cover. Large spatial variations are associated with the 
observations. These are due tO' significant drifting O'f snO'W 
that is characteristic of this basin. The averages of the 
maximum snowcourse and the maximum point measurements 
compared more favourably with the simulations. The syn- 
thesis of the streamflow from the excess snowmelt (figure 
6.2), for the most part, appears to be satisfactory. A 
favourable comparison of observed and simulated flows is 
indicated for the major spring runoff period. However, 
there are appreciable differences between magnitude of flows 
for' the observed and the eimulated, during the periods of 
freeze-up and predominant ice condition of the pre-spring 
thaw period. 

A common set of climatic input data was established for the 

Twenty Mile Creek and Welland River basins. These include a 






five -station average of precipitation and temperature 
inputs. Indices of the other required meteorological 
variables: humidity, net all-wave radiation, estimates of 
potential evaporation and daily wind speed were extrapolated 
from two main meteorological stations (Hornby, (Oakville) 
and Hamilton) that are located nearby in the region. Sub- 
sequently, simulated excess snowmelts serve as the common 
net water input for streamflow simulation in these basins. 

Figure 7.1 shows the simulated indices of snowpack depth and 
density. There are no local snowcourse observations within 
the basins to check on the adequacy of these simulations. 
However, snowcourse observations in nearby basins serve as 
references as to the likely order of magnitude of the simu- 
lated snow accumulation and melt. The simulated depths 
(fig. 7.1) show reasonably good comparison in magnitude and 
trend with those observed at the reference stations. The 
observed reference densities, however, tend to be lower than 
those simulated. Nonetheless, in view of the range of the 
basins' elevation, their nearness to the lakes and the 
inherent large sampling errors that are associated with the 
observations, the simulated values appear to be quite 
reasonable. Figure 7.2 shows the simulated excess melt and 
as an example, the simulated vs observed flows for the 
Twenty Mile Creek basin. A sequence of melt or several 
freeze-thaw cycles, prior to the major spring runoff period, 
are evident. These conditions resulted in difficulties for 
the hydrologic model calibration in matching observed with 
simulated flows. The simulated hydrograph (figure 7.2), 
resulting from the partially calibrated watershed model, 
demonstrates this problem. On the whole, the model is very 
difficult to calibrate in this basin because of the unusual 
physical features of the basin (relatively flat, gentle 
sloping, long and narrow) , The excess melts of the pre- 
spring thaw period contributes primarily to ice build-up in 
the basin and streams. As a result, the total runoff of the 
spring excess snowmelt are substantially increased by the 
water releases from the ice build-up. This is one possible 
explanation for the difficulties in matching observed and 
simulated flows. it is anticipated that further improve- 
ments can be achieved in the application and calibration of 
the hydrologic model for a better prediction of flows in 
relation to the freeze-thaw and ice condition occurrences. 
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Table 2: Monthly Smninary of Historical Daily Streainflow Simulation, 

Wilmot Creek Basin 
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The annual statistics are based on simulation for water-year (Oct-Sept) 
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Figure 5.1: Simulated vs Observed Snowpaek Depth and Density - Wilmot 
Creek Basin: 1969-70 
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Output Sensitivity 

If the simulated excess snowmelts are to be used as 
indices of the net water input to the basin for streamflow 
synthesis, then the extent to which sampling errors in 
precipitation input affect the results of the simulation 
should be a major concern. Simulation experiments were 
performed for a given set of auxiliary meteorological 
variables with differing estimates of basin average pre- 
cipitation input. These include: 2-station average (extreme 
location of gauges) , 5-station average (gauge location 
within the basin boundary), 9-station average (sites, in- 
cluding gauge location outside of the basin) and a 9-station 
average ± 2 standard error of average, respectively. The 
latter two sample groups were considered on the presumption 
that poor gauge (s) or site(s) locations could result in 
either a consistent under or overestimation of the daily 
precipitation input. The 9-station average is assumed to be 
the best estimate of these inputs in the Wilmot Creek basin. 

Figure 8.1, 8.2 and 8.3 display the comparative simu- 
lation results for depth, density and excess melt for the 
winter period of 1969-70. It is seen that the results of 
the simulation for the 5 and 9-station average are not 
significantly different. The result of the 2-station 
average are, similarly, not significantly different, except 
for an approximate 20% less in the maximum depth accumula- 
tion of the pack. On the other hand, the results of the 9- 
station average ± 2 standard error, display dramatic and 
significant differences in magnitude about the mean. The 
underestimation of the precipitation input resulted in about 
a 50% reduction in the accumulated depth and excess melt and 
a shortening of the melt period by about 15 days. The 
overestimation, however, resulted in an approximate 100% 
increase in the depth accumulation, about a threefold 
increase in accumulated excess melt and an extension of the 
melt period by about 21 days. The indications are that a 
reasonable number of carefully selected precipitation sites 
could provide satisfactory estimates of precipitation input 
for the model. Poorly located gauges or sites and erroneous 
precipitation input limit the usefulness of the model in 
simulating snow accumulation and excess melt, for use in 
synthetizing basin runoff. 
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Peak Plow Predictiori 



For the purpose of demonstrating the applicability of the 
snowmelt model in simulating critical melt rates for syn- 
thetizing possible flood flow peaks, hypothetical critical 
sequences of meteorological variables are superimposed at 
maximum snowpack water equivalent storage to create a 
critical melt period. 

A frequency of a 7~day and a 14-day moving average of 
maximum value for each meteorological variable (excluding 
precipitation) , is applied, in turn, to the maximum snowpack 
storage of the 1969-70 snow season. Figure 9 shows the 
predicted flood flows, for the imposed critical melt se- 
quences, as compared to the observed and simulated flows ae 
occurred under normal melt conditions. 



The assumption is that such a procedure provides a suitable 
method for estimating magnitude and distribution of expected 
snowmelt runoff events according to likely encountered 
snowpack environmental conditions of projected critical 
frequency. It is expected that the peak flow and the dis- 
tribution of flood flows will vary according to the choice 
of the frequency or critical sequence of meteorological 
input upon a given snowpack storage condition. However, the 
volume of runoff will remain the same for a given snowpack 
maximum water equivalent storage and rainfall input. 
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SOMMARY AND CONCLUSIONS 

The MOEHYDR- snow accumulating and melt model was 
formulated and tested, in complement with a total watershed 
hydrologic model, for its suitability in augmenting reliable 
simulations of winter flows in Ontario drainage basins. 
Satisfactory simulations of daily variation in snowpack 
depth, water equivalent or density, are achieved for several 
winter periods. Adequate matching of simulated excess 
snowmelt and basin runoff responses are established. As a 
result, satisfactory watershed model calibration and the 
reproduction of historical streamflow are achieved. Problems 
of intermittent freeze-thaw cycles and ice condition build- 
ups that are characteristic experiences in these parts of 
Ontario, created difficulties in the model calibration. The 
snowmelt model output is shown to be very sensitive to poor 
or erroneous precipitation input. However, a precipitation 
sample size in the ord^r of 5 to 9 sites provides adequate 
estimates of the daily precipitation inputs. 

The applicability of the snowmelt model to aid in 
estimating daily flood flows was demonstrated by imposing a 
presumed critical sequence of meteorological variables 
favourable for melting, on maximum snowpack water equivalent 
storage conditions, just prior to the spring runoff period- 
This procedure implies that for any given snowpack at 
maximum storage, different critical sequences of meteoro- 
logical variables may be imposed on the pack to simulate 
varying melt rates and flood flow sequences, thereby imply- 
ing that the snowmelt model can Eerve as a useful design 
tool . 
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